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ABSTRACT 

The X-ray emission from the super-massive star 77 Car is simulated using a three di- 
mensional model of the wind-wind collision. In the model the intrinsic X-ray emission 
is spatially extended and energy dependent. Absorption due to the unshocked stellar 
winds and the cooled postshock material from the primary LBV star is calculated 
as the intrinsic emission is ray-traced along multiple sightlines through the 3D spiral 
structure of the circumstellar environment. The observable emission is then compared 
to available X-ray data, including the lightcurve observed by the Rossi X-ray Timing 
Explorer (RXTE) and spectra observed by XMM- Newton. The orientation and eccen- 
tricity of the orbit are explored, as are the wind parameters of the stars and the nature 
and physics of their close approach. Our modelling supports a viewing angle with an 
inclination of ~ 42°, consistent with the polar axis of the Homunculus nebula (jSmithl 
l2006t ). and the projection of the observer's line-of-sight onto the orbital plane has an 
angle of ~ — 30° in the prograde direction on the apastron side of the semi-major 
axis. 

However, there are significant discrepancies between the observed and model 
lightcurves and spectra through the X-ray minimum. In particular, the hard flux in 
our synthetic spectra is an order of magnitude greater than observed. This suggests 
that the hard X-ray emission near the apex of the wind-wind collision region (WCR) 
'switches off' from periastron until 2 months afterwards. Further calculations reveal 
that radiative inhibition significantly reduces the preshock velocity of the companion 
wind. As a consequence the hard X-ray emission is quenched, but it is unclear whether 
the long duration of the minimum is due solely to this mechanism alone. For instance, 
it is possible that the collapse of the WCR onto the surface of the companion star, 
which would be aided by significant inhibition of the companion wind, could cause 
an extended minimum as the companion wind struggles to re-establish itself as the 
stars recede. For orbital eccentricities, e < 0.95, radiative braking prevents a wind col- 
Hsion with the companion star's surface. Models incorporating a collapse/disruption 
of the WCR and/or reduced preshock companion wind velocities bring the predicted 
emission and the observations into much better agreement. 

Key words: hydrodynamics - stars:early-type - X-rays:stars - stars:binaries - 
stars:winds, outflows - stars:individual(77 Carinae) 



1 INTRODUCTION 

The super-massive star 77 Car provides a unique labo- 
ratory for the study of massive star evolution. Situ- 
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ate d in the Carina nebula, at a distance of 2.3 kpc 
(Davidson fc Humphrevj|l997l ) , rj Car has been identified as 
a luminous blue variable star (LBV), a short-lived (~ 10*yr) 
phase in the life of a massive star occuring after the main 
sequence and preceding the Wolf-Rayet phase. 77 Car is be- 
lieved to have a current mass of ~ 80 — 120 M© , and an ini- 
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tial mass > 150 Mq (iHillier et al.ll2001^ . Over the past 200 
years episodes of rapid mass loss have been observed from 
the star, and in 1843 the "Great Erupti on" ejected more 
than 10 Mq of matter jSmith et al.ll200^ from the central 
object over a ~ 20 year period, forming the bipolar nebula 
known as the Homunculus. In the 1890's a further smaller 
eruption w ith a mass loss of 1 Mq produced the Little Ho- 
munculus jlshibashi et al.ll2003h . The central object is now 
shrouded by a large amount of dense absorbing gas and hot 
and cold dust which makes observations at optical and UV 
wavelengths difficult. 

Although many questions still remain about the na- 
ture of this system, a periodicity of ~ 5.5 yr in the X- 
ray, radio, mm, IR a nd optical is now clearly established 
(jPamineli et al.ll2008l . and references therein), and has pro- 
vided strong evidence for the presence o f a binary compan - 
ion. With an orbital eccentricity of ~ 0.9 (jSmith et al.ll2004 ). 
the stellar separation at periastron may only be ~ 1 .6 au 
(with a semi-major axis of 16.64 au, Imilier et al]|200ll l. The 
radius o f the primary s tar is unclear but it may be as high 
as 1 au ()Daminelilll996h . The companion star has remained 
very elusive, but is thought to have a mass ^ 30 Mq , and 
to be either an O star or a WR star ijPittard fc CorcoranI 
|2002| . hereafter PC02). In this paper the more massive star, 
the LBV, will be referred to as the primary, and the smaller 
companion will be referred to as the secondary. Typical or- 
bital and stellar parameters for the binary system are noted 
in Table [U 

While the binary model is now commonly accepted, 
the orientation of the system, the mass-loss rates of 
the stars, and the key physics occuring around perias- 
tron passage are much more controversial. Most works 
favour the secondary star moving behind the primary 
around pe riastron (e.g. 'Daminoli"l99fr; 'Pittar d et al.l 19981 : 
Corcoran ct al, ^01; Corcoran .200.5; Akashi et al.l l2006l: 
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20081: lOkazaki et all 



20071: [Nielsen eV al.1 l2007i : iHenlev et all 



20081 ) ■ though some others favour a 



system orientated such that the secondary star is posi- 
tioned in front of the primary during pe riastron passage 
(|Falceta-Gori(;aIves et aLll2005l: IXashi fc Sokcr.. 20071 ). or at 
quadrature ( Ishibashl 200ll : Smith et al.l 2004l l . Estimates 
of the mass-loss rate of the primary star range from a 
few times 10"* (|Corcoran et al.l 1200 ll : IPittard fc CorcoranI 

" 200ll : 



20021) to a few times 1 0"^ M, 



©yr 



( Hillier et al.l 



van Boekel fc et~aLl |2003| ). Observations of the reflected 



emission from the Homunculus by ISmith et al.l (|2003l ) re- 
vealed a higher degree of absorption at high latitudes, im- 
plying a higher rate of mass-loss from the primary star near 
the poles. In contrast, the X-ray emission from the WCR 
is most sensitive to the mass loss near the orbital plane. If 
the orbital plane is aligned with the equatorial region of the 
winds this may explain some of the difference in the esti- 
mated mass-loss rates. 

The X-ray emission has proved to be extremely use- 
ful for deciphering the complex phenomena of rj Car, be- 
cause it is relatively un-attenuated and can be used to 
probe deep into the centre of the system. The keV X-ray 
emission displays a dominantly thermal spectrum, and in- 
dicates the presence of hot shocked gas at a temperature of 
~ lO'^ - lO'^K. While there are several possible methods for 
generating high temperature plasma, the X-ray lightcurve 
provides strong support to the binary scenario, where the 



X-ray emission naturally arises from the high- speed collision 
of the winds in a massive binary system (e .g. Istevens et al.l 
ll992l : |Pittard fc Stevens^ll997l : IPittardll2007^ . Since the LBV 
wind is slow, and cannot be shock heated to high tempera- 
tures, the X-ray emission in this scenario must arise from a 
much faster companion wind. 

Fits to a grating spectrum from the Chandra satel- 
lite obtained mass loss rates and terminal wind speeds of 
2.5 X lO~*M0yr~^ and 500kms~^ for the primary star, 
and 1.0 X lO'^M^yr"^ and 3 000 km s~^ for the secondary 
star (|Pittard fc Corcoran|[20o3 ) . The spectrum was attained 
near apastron, when 2D hydrodynamic simulations (which 
ignore orbital motion) give a good description of the WCR. 

A particularly interesting feature of the X-ray emission 
is the long and deep minimum which coincides with perias- 
tron passage of the stars. During this interval emission from 
the system undergoes dramatic variations across a broad 
range of wavelengths, with high excitation lines displaying 
a rap id reduction or "collapse" in intensity (|Damineli et al.l 
|2008| ). We construct a 3D model of the WCR in which rapid 
and large changes in the structure occur during periastron 
passage. The model allows us to investigate the X-ray emis- 
sion throughout the entire orbit, so that the orbital orienta- 
tion, the eccentricity, and the wind parameters of the stars 
can be determined. We then explore the physics of the WCR 
during periastron passage. The layout of this paper is as fol- 
lows: § [2] discusses a number of scenarios which have been 
suggested to explain the behaviour of the X-ray minimum; 
§ [3] contai ns a brief description of the dynamic model de- 
veloped by iParkin fc Pittard! (|2008l . hereafter PP08) which 
we apply to r] Car in § |4l § [5] examines a number of mecha- 
nisms which may be important for bringing about the X-ray 
minimum; and § [6] concludes and outlines possible future 
directions. 



2 THE X-RAY MINIMUM 



2.1 Observations 



The Rossi X-ray Timing Explorer (RXTE) has observed 
77 Car since 1996, so that data over more than 2 orbital 
cycles is now available. The X-ray lightcurve is shown in 
Fig. [1] Following apastron, there is a gradual increase in 
emission over approximately 2 years, followed by a rapid de- 
cline to a deep minimum which lasts for approximately 60 
days (~ 0.03 in orbital phase), then a shallower egress out 
of the minimum to a roughly constan t luminosity prior to 
its gradual rise again in the next cycl e (|lshibashi et al.|[l999l : 
ICorcoran et"alll200ll : [Corcoranll2005l ). 

X MM-Newton spectra analysed by iHamaguchi et al] 
(|2007l ) showed that the flux in the 2-10 keV range at the 
start of the minimum was 0.7% of the maximum value ob- 
served by RXTE before the minimum. Suprisingly, the flux 
in the latter half of the minimum increases by a factor of 
5 from the lowest observed value, indicating that the mini- 
mum has two states. Some of the observed X-ray flux during 
the minimum is emission from an earlier orbital phase w hich 
is scattered by the Homunculus (jCorcoran et al.ll2004h . 
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Figure 1. The 2-1 keV X-ray lig htcurvc of rjCar corrected for instrumental background over two pcriastron passages taken with 
the RXTE satellite l lCorcorarJl2005l) . The data from ~ 1996 — 2003.5 (dotted) has been super-imposed over the 2002-present data for 
comparison. For the majority of the orbit the flux is broadly constant. Prior to pcriastron passage there is a distinctive rise to maximum, 
then a sharp fall to a deep minimum. There is a noticeable difference between the peak magnitudes prior to the minimum of the last 
two cycles. The dotted horizontal line is an estimate of the flux due to cosmic background which contaminates the RXTE/PCU2 field 
of view. 



2.2 Possible explanations 

There are currently a number of possible explanations for 
the minimum in the X-ray lightcurve which may, or may 
not, be mutually exclusive: 

• The eclipse model: The region of the WCR responsible 
for emitting 0.1-10 keV X-rays may be fully occulted by the 
primary star itself, or by its dense, optically thick wind (i.e. a 
wind eclipse) if the secondary moves behind the primary star 
at periastron. The lack of variation in the plasma tempera- 
ture measured in the XMM-Newton spectra during the min- 
imum is supportive of an eclipse of the source, in which case 
residual emiss i on ar ises from regions further downstream. 
lOkazaki et al.l (|2008l ) find that pile-up of the dense primary 
wind after periastron can cause an extended minimum via a 
wind-eclipse. On the other hand, the deformation of the Fe 
XXV profile and the relatively weak fiuorescence Fe line in- 
tensity during the mini mum may suggest an int rinsic fading 
of the X-ray emissivity (|Hamaguchi et al]|2007h . 

• Increased mass-loss: LBV's exhibit variab ility on a 
number of timescales (e.g. Ide Groot et ahlbooj ), and are 
known to go through eruptive phases where large amounts 
of stellar material are ejected into the surrounding in- 
terstellar medium. The p rimary star may be undergo - 
ing S Dor- like fluctuations (|van Genderen fc Sterkenll2007h . 
which may reach a maximum around periastron with 
a shell-ejection event (such an event was proposed by 
iDavidson fc Humphrevslll997l . as an alternative to the bi- 
nary eclipse model for explanation of the X-ray minimum). 

Theoretical work on tidal interactions in binary sys- 
tems shows that tidal oscillations m ay result in a phase 
of increased non-sp herical mass- loss ( Moreno et al.l l2005l: 
iKashi fc Sokeij[2008l . and references therein) . Corcoran et al.l 
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Figure 2. The large scale structure of the contact discontinuity 
(which separates the two winds) viewed from above the orbital 
plane. The white region marks the unshocked wind of the primary 
star. In this model the unshocked secondary wind is entirely con- 
tained within the coloured region. The primary star is located at 
the origin and the companion star orbits about this point. The 
WCR forms a spiral structure due to the orbital motion of the 
stars. Asymmetry in the distribution of the unshocked compan- 
ion wind is caused by the highly eccentric (e = 0.9) orbit. The 
thinning of the coloured region towards the tail of the spiral is an 
artifact of the model (see PP08), due to the fact that the WCR 
is traced back for only two orbits. The parameters used in this 
simulation are noted in Table [T] 



4 E. R. Parkin, J. M. Pittard, M. F. Corcoran, K. Hamaguchi, and I. R. Stevens 



Table 1. Assumed parameter s for r) Car, r) = M2Vao2l ^li^oa l is t he wind momentum ratio and a is the semi-major a xis of the orbit. 
Refere nces are as follows: 1 = IPavidson fc Humphrevsi l ll997t) . 2 = iHillier et al] 1 I2OOII) . 3 = iPittard fc CorcoranI l[20o2), 4 = ICorcoranI 
l|2005l '). 5 = IOorcoran fc Hamaguchil l|2007t )." After an extensive parameter space exploration wo found that a composite (ISM + nebula) 
column density (between the observer at Earth and the edge of the simulation box, 1500 au from the stars) of 1 X 10^^ cm~^ provided 
better fits to the X-ray flux in the s oft band (2-5 keV). This value is consistent with, but slightly lower than, the value of 5 X lO'^'^ cm^^ 
derived bv lHamaguchi et al.l |20o3) for the "the absorption beyond the central source", which likely represents the absorption column 
between the secondary star at apastron and the observer at earth, ij* is taken to bo the radius of the gravitationally bound region of 
the star; for the primary star the photosphere will exist somewhere in the stellar wind. 
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Figure 3. The structure and position of the WCR around periastron passage. The colour scale shows the surface density of the postshock 
primary wind in gcm~^. The primary star is shown as a green circle. The tick marks on the axis mark a distance of 5 au (note the 
different scale of each panel). From left to right the plots cor respond to 20 days b efore periastron, periastron itself, and 40 and 81 days 
after periastron, and compare well against the SPH model of lOkazaki et all bOOSi ). 



(|200ll ) found that incorporating a factor of ~ 20 increase in 
the primary star mass-loss rate for 5<j) ~ 0.04 improved the 
fit to the duration of the minimum. High velocit y absorption 
lines ( ~ 750 kms"^ for He I A6678) observed bv lStahl et al] 
(|2005l ) around periastron in the reflected emission from the 
Homunculus nebula may be evidence of a period of increased 
mass- loss. 

• Shut-down of the companion's wind: If the WCR moves 
into the acceleration region of the secondary's wind (due 
to contraction of the stellar separation), a stable momen- 
tum balance may be lost, resulting in the primary's wind 
overwhelming the secondary's and causing the WCR to col- 
lapse onto the surface of the seco ndary star. This may be 
aided by th e radiative inhibition (IStevens fc Pollocklll994l ') 
or braking (|Gavlev. Owocki fc Cranmed I1997D " of the sec- 
ondary star's wind at periastron due to the high stellar lumi- 
nosity of the primary star (~ a few xlO^ L0). Evidence for 
a collapse of the WCR onto the secondary star exists in the 
behaviour of the high ly ionised lines around the minimum 
JPamineli et al.|[2008h . This possibility, and the potential oc- 



curance of accretion of the primary wind by the companion 
star, was explored bv lAkashi et al.l (|2006l ). 



3 THE DYNAMIC MODEL 

To successfully model the X-ray minimum the motion of 
the shocked and unshocked gas relative to the stars must 
be taken into account. Unfortunately, with current compu- 
tational resources, it is difficult to explore a wide range of 
parameter space with 3D models, even when using sophisti- 
cated grid or particle based hydrodynamic codes. The high 
orbital eccentricity of 77 Car compounds this problem, and 
3D simulations of 77 Car have only examined a small region 
of pa rameter space to date (|Pittardl l200d : lOkazaki et al.l 
|200q). While both of these efforts revealed important in- 
sights into the dynamics of the WCR, only a small part 
of the o rbit a round periastron passage was simulated in 
IPittardI (|2000l ). while the assumption of isothermality in 
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lOkazaki et alj (|2008l ) does not allow direct calculation of 
the X-ray emission. 

In this paper we use a recently developed 3D dynamical 
model which allows efficient exploration of parameter space. 
A detailed description can be found in PP08 where a demon- 
stration with lower eccentricity O-l-O and W R-|-0-star bi- 
naries is given. The model uses equations from lCanto et all 
1 19961 ) to determine the ram pressure balance and position 
of the contact discontinuity (CD) between the winds (modi- 
fied by an aberration angle due to orbital motion). The flow 
downstream of the "apex" of the WCR is then assumed to 
travel ballistically. The motion of the stars causes the con- 
tact surface to wind up and form large scale spiral structures 

i Fig. [H, reminisce nt of the beautiful "pinwheel nebulae" 
TuthiU et aDl2008l '). The shape of the WCR around perias- 
tron passage is highly distorted due to the rapid transit of 
the stars and the slow speed of the primary wind (Fig. [Sjl . 

The X-ray emission from the WCR is a function of the 
gas temperature and density. Since the dynamical model 
does not contain such information we use a grid-based, 2D 
hydrodynamical calculation of an axis-symmetric WCR to 
obtain this. The resulting emission as a function of off-axis 
distance is then mapped on to the coordinate positions in 
the 3D dynamical model. In this way we obtain the benefit 
of effectively modelling the thermodynamic and hydrody- 
namic behaviour responsible for the production of the X-ray 
emission, while simultaneously accounting for the effect of 
the motion of the stars on the large-scale structure of the 
WCR and the subsequent wind attenuation. Since the hydro- 
dynamic calculation is 2D, the computational requirements 
remain low. 

The X-ray luminosity per unit volume is given by 
T{E) = n^A{E,T), where n is the gas number density 
(cm~^) and A{E,T) is the emissivity (erg cm^s"'^) as a 
function of energy E and temperature T for optically thin 
gas in coUisional ionization equilibrium. A{E, T) is obtained 
from look-up tables calc ulated from the MEKAL plasma 
code (|Liedahl et al.|[T995l , and references therein), contain- 
ing 200 logarithmically spaced energy bins in the range 0.1- 
10.0 keV, and 101 logarithmically spaced temperatures from 
lO'' to lO^K. Solar abundances are assumed. 

The orbit of the stars is modelled in the xy plane, so 
that viewing angles into the system can be described by the 
inclination angle that the line-of-sight makes with the z axis, 
i, and the angle the projected line-of-sight makes with the 
major axis of the orbit, 6 (positive values correspond to an 
angle subtended against the positive x axis in the prograde 
direction). The components of the unit vector along the line- 
of-sight, u, are 

Ux = cos 9 sin i, 
Uy = sin sin i, 
= cosi. 

To determine the observed attenuated emission ray-tracing 
is performed through the spiral structure. Specific details of 
the emission and absorption calculations are now discussed. 

3.1 Emission from the shocked secondary wind 

The importance of cooling i n the WCR can be quantified 
using the cooling parameter (jStevens et al.|[T992l ) 



_ tcool _ 'ygrfl2 ^ 

where wg is the preshock wind velocity in units of 
1000 km s~ , di_2 is the distance from the star to the con- 
tact discontinuity in units of lO^^cm, M_7 is the mass-loss 
rate of the star in units of 10~^ yr~^ , tcooi is the cooling 
time, and tcsc is the escape time of gas out of the system 
(~ di2/cs, where Cs is the postshock sound speed), xi and 
X2 are the cooling parameter for the postshock primary and 
secondary winds respectively. 

In the case of rj Car, the postshock primary wind is 
strongly radiative {x — 10~* — 10~^), rapidly cools, and col- 
lapses to form a thin dense sheet. As previously noted, this 
wind is not expected to produce thermal X-rays above 1-2 
keV as the relatively low wind speed {vi ~ 500kms~^) does 
not shock- heat gas to high enough temperatures. In contrast, 
the high velocity of the secondary's wind (w2 ~ 3000 km s~^) 
causes its postshock flow to be adiabatic (x ^ 10) around 
most of the orbit. To first order the entirety of the observ- 
able X-ray emission originates from the postshock secondary 
wind although mixing at the CD due to dynamical insta- 
bilities (PC02; see also Fig. I17|l is a complicating factor. 
Note that lower values of X2 at periastron indicate that 
cooling may be important for a short period (Fig. I16p . es- 
pecially if the secondary wind is clumpy and/or if radia- 
tive inhibition reduces its preshock speed. These possibili- 
ties are investigated in more detail in §[S1 In the meantime, 
we assume that the shocked secondary's wind is adiabatic 
throughout the entire orbit, emits the entirety of the X-ray 
emission, and that the intrinsic X-ray lumin osity scales as 
l/dsep, where daep is the stel lar separation (jStevens et al.l 
1 19921 : iPittard fc Stevenslll997l ) . 

3.2 Attenuation 

The intrinsic X-ray emission is attenuated by the unshocked 
winds. In wide binaries one can usually assume that the 
winds are instantaneously accelerated. However, in rj Car the 
primary's wind may accelerate very slowly and this may 
affect the degree of attenuation observed. When describing 
the wind acceleration using a /3-velocity law, a typical value 
fo r the slow acceleration of an LBV wind would he P — 
4 (IBarlow fc Cohen|[l977l : iPauldrach fc Pulslll990l ). To test 
whether this is an important factor we compared simulations 
adopting accelerating and instantaneously accelerated winds 
and found that the difference in attenuation was small at all 
phases, and for all lines of sight. 

Attenuation by the shocked winds also occurs. Since 
the shocked secondary wind is assumed to be hot, it is also 
assumed to have negligible attenuatiorQ. In contrast, the 
dense, cool layer of postshock primary wind is a signifi- 
cant source of absorption. The surface density of this layer 
changes by over two orders of magnitude between periastron 
and apastron (a oc d~^), resulting in a dramatic increase in 
attenuation around periastron for X-rays. For viewing angles 
closely tangential with the shock surface the combination of 

^ This is a suitable assumption for the purposes of our model, 
however it has been suggested bv lHenlev et al.l bOOSi ) that atten- 
uation by the companion wind plays a role in shaping the emission 
line profiles. 
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Figure 4. The effect of various attenuation mechanisms on the observed 2-10 keV emission as a function of orbital phase over an entire 
orbital period (left), and over the orbital phase range 0.90-1.10 (right). The un-attenuated intrinsic emission (I) is shown by the solid red 
curve. Various attenuation mechanisms arc then sequentially added: occultation (O), wind attenuation (W), and attenuation through 
the thin dense postshock primary wind layer (SD). The viewing angles for these curves are i = 90° and 6 = 0°, with the other model 
parameters listed in Table[T] The aberration of the WCR due to orbital motion and interstellar -|- nebula absorption (1 X lO^^cm^'^) are 
included. 
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Figure 5. A 3D plot of the occultation of the WCR by the pri- 
mary star at periastron {ip = 0.0) with e = 0.9 and a viewing 
angle of i = 90° and 6 = 0°. The region occulted by the primary 
star (displayed in black) is shown in green. Although the primary 
star is directly in front of the secondary star for the chosen values 
of i and 9, due to aberration the shadow does not fall directly on 
the apex of the WCR. 



the increased postshock gas density and the pathlength de- 
pendent absorption described in PP08 is sufficient to fully 
absorb X-rays with energies < 4 keV. 

3.3 Observed emission 

Prior to modelling the RXTE lightcurve, we explore the 
effect on the observed X-ray lightcurve of including various 
occultation and attenuation mechanisms. A lightcurve from 
a model with the parameters listed in Table [1] and viewing 
angles i = 90° , S = 0° is shown in Fig. |4l The skew of 
the WCR due to orbital motion causes the minimum to be 
asymmetric, as observed. A particularly interesting feature 



is the shoulder seen on the egress out of minimum. This is 
caused by attenuation through the dense layer of shocked 
primary wind (see below). 

Occultation by the primary star (Fig. [5]) causes a consid- 
erable reduction in emission around periastron though the 
effect is very short-lived and does not cause the observed 
minimum (see Fig. |4]). Instead, the largest source of atten- 
uation arises from the unshocked stellar winds except for a 
phase interval of ~ 0.01 during the minimum when absorp- 
tion through the shocked primary wind dominates (Fig. |6]). 
The smallest values of the emission weighted column (PP08) 
occur at apastron, when the sight lines from the head of 
the WCR initially pass through the secondary wind. A cir- 
cumstellar column of 2 x 10^^cm~^ (the absorption from 
shocked and unshocked winds within the simulation box, 
see PP08), plus a column of 1 x lO^^cm"^ which includes 
absorption in the Homunculus and ISM, is sufficient to ab- 
sorb ~ 99.99%, 90%, and 20% of the 1, 2, and 5 keV X-rays 
from the WCR respectively. The column density begins to 
increase as the stars approach periastron, then there is a 
dramatic rise in the column through the unshocked winds 
at (/) ~ 0.99 which highlights the shift of the lines-of-sight 
into the high density primary wind. This is accompanied 
by a peak in absorption by the shocked primary wind, rep- 
resented by the surface density (SD) component. The two 
peaks in the SD component column at phases cj> ~ 0.99 and 
1.02 are produced by sight lines originating from points with 
the highest intrinsic luminosity becoming closely tangential 
with the shock surface. The short duration in phase of the 
peak in column density (which briefly reaches ~ lO'^^cm"'^) 
denotes the rapid motion of the stars through periastron 
passage. 

The attenuation along lines-of-sight close to the positive 
X-axis (i.e. i = 90°, d = 0°) at phases near apastron is 
composed of large volumes of unshocked secondary wind and 
slices of dense primary wind (c.f. Fig. [2]), with the former 
dominating. 

Increasing the radius of the primary star causes the 
depth, and width, of the occultation component of the min- 
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Figure 6. The effect of various attenuation mechanisms on the emission weighted column density as a function of orbital phase. The 
total column (red), attenuation due to the unshocked stellar winds (blue), and attenuation due to X-rays intersecting the dense, cold 
layer of postshock primary wind (green) are shown over an entire orbital period (left) and over the orbital phase range 0.90-1.10 (right). 
Interstellar -|- nebula absorption (~ 1 X 10^^ cm~^) has not been added to these plots. The viewing angles are i = 90° and 9 = 0°. 
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Figure 7. Synthetic 2-10 keV lightcurves for various inclination angles, i, over an entire orbital period (left) and over the orbital phase 
range 0.90-1.10 (right). An inclination angle of 45° provides the best resemblence to the observed morphology (Fig. [TJ. The parameters 
used to calculate these results are noted in Table [T] 6 = 0° in each of these models. 




Figure 8. Synthetic 2-10 keV lightcurves for an inclination i = 90°, eccentricity, e = 0.9, and various line-of-sight angles, 6, over an 
entire orbital period (left) and over the orbital phase range 0.90-1.10 (right). All other parameters are noted in Table[T] 
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Table 2. Parameters used for the models shown in Figs. [TT] and 
1121 In all simulations the terminal velocity of the primary, Vool, 
and secondary, Voo2, stars wind's are 500 and 3000 kms"'^ re- 
spectively. 
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imum, and the actual minimum to increase slightly, but the 
high eccentricity combined with the rapid orbital motion 
of the stars around periastron results in very little change 
to the attenuating column densities. The lightcurve is also 
relatively insensitive to the rate of the acceleration of the 
primary wind (which affects the density near to the star) 
and to the velocity cut at which the flow in the WCR is 
assumed to behave ballistically (see PP08). 



4 RESULTS 

4.1 Varying model parameters 

In the following subsections we explore the effect of varying 
the orientation of the orbit, the wind momentum ratio, and 
the eccentricity. 

4.1.1 Variation with orbital inclination 

Fig. [7] shows how varying the orbital inclination angle 
changes the observed lightcurve. For values of i = 60° — 90° 
the synthetic lightcurves are largely similar. For i ~ 45°, 
similar to the 42° inclination angle deri ved for the polar 
axis of the Homunculus nebula (Smith 20061 ). the ratio of 
the pre- and post-minimum X-ray emission is reduced con- 
siderably and is in much better agreem e nt with the RXTE 
data (see also ISteiner fc Damineiill2004l: lAkashi et al.l l 2006l : 
iKashi fc Sokeill2007l : lOkazaki et al.ll2008l '). Reducing the in- 
clination angle further to i = 30° causes the pre-/post- 
minimum flux ratio to become too small. Although reducing 
i decreases the direct occultation of the region of the WCR 
with the largest intrinsic luminosity around periastron pas- 
sage (Fig.O, the wind absorption increases as a larger frac- 
tion of the X-rays initially pass through the primary wind 
(cf. Fig.0. 

4. 1.2 Variation with the orientation of the orbital 
semi-major axis 

As mentioned in § [1] there are competing suggestions for the 
orientation of the semi-major axis. Fig. [8] shows the effect 
on the observed X-ray emission of varying the line-of-sight 
angle, 9, while keeping the inclination angle constant at i = 
90°. In each case, the maximum emission occurs while the 
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Figure 9. Synthetic 2-10 keV lightcurves over the orbital phase 
range 0.90-1.10 for rj = 0.24, an inclination i = 90°, line-of-sight 
angle 9 = 0°, and varying values for e. The orbital and wind 
parameters used are noted in Table [T] 
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Figure 10. Synthetic 2-10 keV lightcurves for ri = 0.24, line-of- 
sight angle 8 = 90°, eccentricity e = 0.9, and varying inclination 
angle, i, over the orbital phase range 0.90-1.10. The orbital and 
wind parameters used to calculate these results are noted in Ta- 
ble [T] Comparison with Fig. [T] shows that these curves do not 
reproduce the observed lightcurve. 



sight line to the apex of the WCR is through the less dense 
secondary wind. 

4.1.3 Variation with orbital eccentricity 

For the majority of the calculations an orbital eccentricity of 
e = 0.9 was adopted. The eccentricity affects the orbital ve- 
locity of the companion star around periastron, the distance 
of closest approach of the stars, and the amount of time the 
companion star and the region of the WCR with the high- 
est intrinsic luminosity (assuming Lx-^„^ ~ but see §[5] 
for caveats) spends immersed in the dense primary wind. 
Lowering the value of e increases the width of the minimum 
towards cj> ~ 0.03 which is desired (Fig.|9}. However, the de- 
cline to minimum then occurs at earlier phases. This can be 
countered by increasing the value of 6 (c.f. Fig.[S|, although 
there are then implications for the pre-/post-minimum lu- 
minosity ratio, and the morphology of the egress out of 
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minimum. On the other hand, increasing the eccentricity 
decreases the duration of the minimum and increases the 
pre-/post-minimum luminosity ratio. Overall, a best-fit is 
obta ined with e ~ 0.9 , in agreement with previous work 
(e.g. iHenlev et al]|2008l 'l. However, if the WCR collapses at 
periastron (see § [5| the eccentricity may be higher. 



4- 1-4 Variation with wind momentum ratio 

Varying the momentum ratio between the winds (by ad- 
justing All) affects the opening angle of the WCR and the 
spatial distribution of the unshocked winds, and changes 
the attenuation and observed X-ray emission, as shown in 
Fig. 111! Decreasing 77 causes the sight lines to the apex of 
the WCR to move into the denser LBV wind at earlier or- 
bital phases, and produces a wider minimum at periastron. 
When r; — 0.024 the opening angle of the WCR is reduced 
so much that the sight line from the WCR apex moves into 
the primary wind at such an early phase that a significant 
rise in X-ray luminosity before periastron is prevented. An- 
other effect of reducing 77 is that more of the secondary wind 
becomes shocked in the WCR. This means that the luminos- 
ity at apastron initially increases as 77 decreases from 0.24. 
However, the opening angle of the WCR soon decreases to 
such an extent that sightlines to the apex of the WCR move 
into the primary wind, so that the luminosity at = 0.5 
declines for 77 < 0.04. 

With Ml = 3.5x IO^^Mq yr"^ (and rj = 0.024) the peak 
column density at periastron increases to ~ 2 x 10^^ cm~^ 
(Fig. I12|l . The X-ray minimum is flat and has approximately 
the correct width {Sep ~ 0.03, corresponding to 60 days), yet, 
the decline to, and egress out of, minimum are too shallow, 
and the rest of the lightcurve is a poor match to the ob- 
servations. If increased mass-loss is responsible for the flat 
bottom of the minimum. Fig. Illl shows that it must be con- 
flned to a short pha se interval around periastron (see also 
ICorcoran et al.|[200ll ). The effect of a short-lived episode of 
increased mass-loss is discussed further in § [S] 



4.2 Comparison to the X-ray data 

In the following subsections we examine flts to the RXTE 
lightcurve and hardness ratio, and XMM-Newton spectra. 



4.2.1 Fits to the RXTE lightcurve 

Having explored how the lightcurves vary with i,0,e and 
77, we performed further simulations to o btain a best flt 
to the RX TE data. The binary models of [Abraham et al.l 
(|2005D and lHenlev et all (|2008l ) suggest i ~ 90° and 70° re- 
spectively. However, for such values we find a strong dis- 
agreement between the RXTE and synthetic lightcurves, 
mainly due to the large difference between the pre- and 
post-minimum X-ray luminosities in the model lightcurves. 
A comparison of the 8 — 180° curve in Fig. |8] and the 
RXTE data shows that an orientation where the secondary 
star is in front of the primary during periastron p assage 
iFalceta-Goncalves et apl2005l: iKashi fc Sok"eill2007h is not 



sensible inclination angle when 8 = 90°, as shown in 
Fig. 1101 Instead, the RXTE data requires the observer 
to be situated close to the semi-major axis with S ~ 0° 
19961: 



200 ll: 



Daminelil 



Pittard et al.l 19981: Corcoran et al.l 



Corcoranll2q05l:lAkashi et al.ll20q6l: iHerilev et al.ll200S 



supported. ISmith et all (|2004h suggested that 8 = 90°. 
However, the required pre-/post-minimum fiux ratio, and 
duration of the minimum, cannot be obtained with any 



Haniaguchi et all l2007l : iNielsen et all l2007l : iDamineh et all 
20081 ) ■ whereby the companion star moves behind the pri- 
mary at periastron or just after. A detailed analysis of pa- 
rameter space using "chi-by-eye" yields be st-fit values of 
i^_42° and 8 ~ 20°, in good agreement with lOkazaki et"al] 
(|2008h . 

A good match to most of the lightcurve is obtained with 
i = 42°, 8 = 0-30°, Tj = 0.18 and e = 0.9 (see Fig. [T3)) . 
though a flat extended m inimu m cannot be prod u ced. I nter- 
estingly, ICorcoranI (120051 ) and iHamaguchi et all (|2007l ) dis- 
cussed two phases to the minimum; an initial deep minimum 
and a shallower phase which begins approximately halfway 
through. This behaviour is seen in the model lightcurves and 
can be understood by examining Figs.[6land ll2l which show 
a step in the emission weighted column density due to ab- 
sorption by the cooled postshock primary wind exceeding 
the absorption through the unshocked winds. Note that this 
step is not visible in the 77 — 0.024 lightcurve in Fig.[TT]as the 
attenuation due to the unshocked winds always dominates. 

The similarity between the model and RXTE 
lightcurves with 77 = 0.18 shows that the system can be 
described well for the majority of the orbit with one set 
of parameters. However, the exact shape of the minimum 
could not be matched with any single set of mass-loss rates 
and terminal wind speeds. In particular, the absorption peak 
(Fig. I12p does not have the necessary duration to create 
a flatfish minimum. The steep rise out of minimum ob- 
served by RXTE also cannot be replicated. This indicates 
a need for some additional physics around periastron pas- 
sage which causes the column to stay at lO^^cm"^ until 
(j) ^ 1.03 {XMM-Newton observations find the column to be 
> 10^^ cm^^ over this phase range) or some physical process 
to reduce the intrinsic emission. A number of mechanisms 
which could effect the structure of the WCR and the ob- 
served emission are discussed in § [S] 

Since the intrinsic X-ray luminosity scales as Lx oc 
M^ in our model (where it is assumed that the post- 
shock secondary wind is adiabatic), the normalization of the 
lightcurve can indicate the mass-loss rates of the stars. We 
find that when model lightcurves produced with the wind 
parameters in Table [T] are compared to the RXTE data the 
normalization of the model is a factor of ~ 2.2 too low, 
which for fixed 77 corresponds to an underestimate in M2 of 
~ 1.5. Various factors contribute to this discrepancy, such 
as the distance to the source and the interstellar column 
adopted in the calculations compared to those of PC02, 
who assumed a distance to 77 Car of 2.1 kpc, and obtained 
a combined circumstellar and interstellar column density of 
~ 7.7 X lO^^cm"^. In contrast, this work adopts an inter- 
stellar -I- nebula column of 1 x lO^^cm"^ (the majority of 
which being due to the Homunculus nebula as the ISM col- 
umn density is only ~ 3 x lO^'^cm"^) and a distance of 2.3 
kpc (determined from long-slit spec troscopic observations 
with the HST bv lffillier et aLlbOOll ). In addition, the lack 
of spatial resolution with RXTE will lead to additional 2- 
10 keV flux from nearby, unresol ved sources, includ i ng th e 
constant components identified bv IHamaguchi iiai]i^Q3). 
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Figure 11. Variation of the attenuated 2-10 keV X-ray emission for e = 0.9, i = 42° , 6 = 20° and varying values of the wind momentum 
ratio rj shown over the whole orbit (left) and over the orbital phase range 0.90-1.10 (right). The orbital and wind parameters are noted 
in Tables [1] and m 
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Figure 12. Variation of the emission weig hted column density (cm-^) for 7? = 0.24, 0.18, 0.12, and 0.024, with e = 
shown over the whole orbit (left) and over the phase range 0.9-1.1 (right). The interstellar -|- nebula column 
additional to these plots. 
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Figure 13. Variation of the attenuated 2-10 keV X-ray emission for rj = 0.18, e = 0.9, i = 42°, and values of 9 of —10° to 30° in 
increments of 10° over the whole orbit (left) and over the orbital phase range 0.90-1.10 (right). The orbital and wind parameters used 
are noted in Tables [T] and [2] The model results are plotted against the second orbit RXTE data (which has a higher pre-minimum flux 
than the first orbit data, see e.g. Figs. [Tl and 1241 for comparison. 
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Figure 14. Variation of hardness ratio with phase for the best fit ?; = 0.18 model and the RXTE data over the 2.5 orbits. The model 
hardness ratio has had the detector response folded in for each specific observation. Due to fluctuations in the calibrated channel energy 
gains there is some jagged structure to model hardness ratio leading up to periastron passage. The large drop at rf> = 1.23 corresponds 
to a gain reduction in the PCU on-board RXTE, and is also incorporated in the folded model. The change at </) = 0.7 is due to a 
change in the detector response model. The soft band (s) is 2-5 keV, the hard band (h) is 7-10 keV, and the ratio is calculated as 
hr = (h — s)/{s + h). The dip in the model ratio at periastron is due to occultation of the WCR by the primary star. Over the orbital 
phase range 0.924 ^ <p ^ 1.023 and 1.924 ^ </> ^ 2.023 the hardness ratio values are from the individual models with lower companion 
wind velocities, discussed in § 15.21 are shown (red). In the 2.0 ^ ip ft 3-0 cycle the model appears softer than in the previous cycle, 
however this is the result of processing the model data with the detector response. 



In order to compare like-with-like a background spectrum 
of these unresolved sources has been added to all of the 
model lightcurves where a direct comparison is made against 
RXTE data, e.g.. Figs. [l3 \M and[24l Despite this addi- 
tion, a small discrepancy remains, and slightly larger M's 
(Ml = 4.7 X 10"* Mq yr"^ and M2 = 1.4 x 10"* Mq yr"^ ) 
than determined by PC02 are required to bring the model 
results into agreement with observations. 

4-2.2 The hardness ratio 

A comparison of the hardness ratio from the best-fit 77 — 0.18 
model and the RXTE data is shown in Fig. ll4l We find there 
is a good agreement between the model and the RXTE data 
for large parts of the orbit. The model hardness ratio ini- 
tially rises towards periastron as stronger absorption pref- 
erentially reduces the soft-band flux. This rise is also seen 
in the RXTE data. However, the most noticeable feature in 
the RXTE data is a sudden softening of the flux which lasts 



throughout the X-ray minimum. In contrast, the model flux 
is hard for the duration of the X-ray minimum, except for a 
brief period when the head of the WCR is occulted by the 
primary star. The lack of agreement between the data and 
the model indicates that the X-ray minimum is not solely 
caused by absorption and occultation effects. Instead, there 
appears to be a substantial intrinsic reduction of the hard- 
est emission in the RXTE data during the minimum. As we 
shall see, this also helps to explain the discrepancy which 
exists between the lightcurves in Fig. 1131 The observational 
data used to calculate the hardness ratio includes emission 
from prominent background sources in the RXTE fleld-of- 
view such as WR 25 and Trumpler 14 and 16. The flux 
from these sources is predominantly in the soft band and 
accounts for the fact that although the = 1.009, 1.015, 
and 1.023 spectra from 77 Car shown in Fig. 1151 clearly have 
higher hard band flux than soft band, the softening in the 
RXTE data may be due to imperfect accounting for the soft 
X-ray background. 
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An interesting result is that there is a softening in the 
observed flux from r] Car after phase ^ 2.1, which is high- 
lighted by the fact that although the RXTE data appears 
very similar through both cycles there is a noticeable differ- 
ence between the model and the data when both have the 
detector response included. This is suggestive of an intrinsic 
variation in the source of the emission, but uncertainties in 
the assumed PCU detector response in the 2-10 keV band 
cannot yet be ruled out. 



4.2.3 The X-ray spectra 

A series of XMM-Newton observations obtained through 
the 2003 perias t ron passage have been analyzed by 
iHamaguchi et ahl (|2007l ). The observed spectra are a com- 
bination of a variable hard component originating from the 
WCR, and a number of other non- variable components. The 
XMM-Newton spectra in Fig. llSl have the non- variable emis- 
sion components from the outer ejecta, the X-ray Homuncu- 
lus nebula, and the central c o nstan t emission component 
identified bv IHamaguchi et al.l (|2007l ) removed. These spec- 
tra are compared to those produced from our best model 
fit to the RXTE lightcurve. The Fe line at ~ 7 keV has 
a large fluorescence component which is not modelled in 
the synthetic spectra; the goal of this work is to repro- 
duce the broad-ban d spectra, and not the lines. Recently, 
iHenlev et al.l (|2008h examined the variability of the Sulfur 
and Silicon emission to study the flow dynamics within the 
WCR. Future work will focus on developing the dynamic 
model to perform individual studies of spectral lines. 

At (j) = 0.924 we find that there is a significant de- 
ficiency in the emission from our 'standard' model (with 
Voo2 = 3000 kms~^) at the softer end of the spectrum 
{E = 1 — 4 keV) . Thus the model spectrum is slightly harder 
than the data, indicating that the postshock temperature 
of the companion's wind is over-estimated. The tempera- 
ture reduction in the postshock companion wind could be 
partly caused by mixing of the cold postshock primary wind 
into the hot postshock secondary wind, which softens the 
spectrum, or radiative inhibition (see § I5.2|l . An alternative 
explanation is that ef ficient particle acceleration causes a 
weaker sub-shock (c.f. iPittard fc Doughertvl[2006l ). Strong 
evidence for particle accele r ation in rj Car has recently been 
presented by iLevder et al.l (|2008l ) , who detected a hard X- 
ray tail (£ ~ 10 — 30keV). Including the effects of efficient 
particle acceleration in the models used in this paper is be- 
yond the scope of the current work. The <j} — 0.988 spec- 
trum occurs at the onset of the minimum. Although the 
stars are close at this point, it is still unlikely that the WCR 
has entered the acceleration region of the companion's wind 
(dsep ^ 3.7 au for a = 16.64 au and e = 0.9 c.f. Fig. EO)) . 

The q uasi-periodic " fiares" observed in the RXTE 
lightcurve (jCorcora nl bOOSl ) are not reproduced in the syn- 
thetic data. These may result from dynamic instabilities in 
the WCR or stellar winds and cannot be reproduced in the 
dynamic model used in this work which is assumed to have a 
stable WCR. The XMM-Newton observations at = 0.988 
and 0.990 occur at the bottom and top of a flare peak re- 
spectively, which accounts for some of the discrepancy seen 
between the data and model spectra. 

The spectra from the 'standard' model at = 
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Figure 16. Variation of the cooling parameter, X2, for the com- 
panion wind in the orbital phase range = 0.9 — 1.1. Results are 
shown for calculations with (green) and without (red) the inclu- 
sion of the relative motion of the stars, which affects the position 
of the ram pressure balance and the preshock wind speeds. Prior 
to periastron the stars approach each other, and the relative wind 
speeds are higher, increasing x compared to the static case. The 
opposite effect occurs after periastron. Both curves reach a min- 
imum at periastron. For simplicity the skew of the WCR due to 
orbital motion is neglected, and the wind is assumed to be per- 
fect ly smooth. Clumping and e fficient particle acceleration (see 
e.g. jpittard fc Doughertvl EoOSi ) will reduce X2 below the level 
shown. The effect on X2 due to radiative inhibition is also shown 
(see §[52}. 



1.009, 1.015, 1.018, and 1.023 also show significant excess 
fiux at high energies. 



5 THE NATURE OF THE X-RAY MINIMUM 

In the previous section model simulations were performed 
which provided a reasonable fit to the lightcurve as a whole 
but did not give a good fit to the minimum. A detailed 
examination of the spectrum and hardness ratio around the 
minimum showed some significant discrepancies between the 
model results and the observations, with the model emission 
generally harder than the data during the minimum. In this 
section possible solutions for resolving this discrepancy are 
discussed. 

A plausible scenario is that around periastron the 
preshock speed of the secondary wind is reduced to the point 
that it does not shock to high enough temperatures to pro- 
duce the hard X-ray emission which is normally seen. This 
change in the preshock velocity may be initiated by the 
WCR moving into the acceleration zone of the secondary 
wind due to the high orbital eccentricity, or to a period 
of enhanced mass-loss from the primary star. Alternatively 
the acceleration of the secondary wind may be inhibited to- 
wards the primary star by its enormous radiation field. In 
addition, the shocked secondary wind may become radia- 
tive around periastron (especially if radiative inhibition is 
a strong effect). If both winds strongly cool, it is possible 
that the WCR disrupts and breaks up into a mass of cold 
dense blobs, each surrounded by oblique shocks as they are 
impacted by the supersonic winds (i.e. the WCR "discom- 
bobulates" - Kris Davidson, personal communication). The 
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Figure 15. 1-10 kcV model spectra with 77 = 0.18, i = 42° and 9 = 20° plotted against XMM-Newton data (see iHamaeuchi et al.l[2007t ) 
as a function of phase. In the 'standard' model the companions wind is assumed to have a preshock speed of 3000 kms~^. Results are 
also shown for models with a reduced preshock companion wind velocity. A model which includes a collapse/disruption of the WCR is 
also included in the panels at phases (p = 1.009, 1.015, 1.018, and 1.023. Models with lower wind speeds or a collapse of the WCR match 
the observed flux above 5 keV during the minimum much better. 
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fast companion wind may then be slowed through a sequence 
of oblique, radiative shocks, each with a lower postshock 
temperature than would be obtained behind a single global 
shock. 

However, the above processes by themselves cannot ac- 
count for the long duration of the minimum, and the real 
picture is likely to be even more complicated. There is a 
possibility that the WCR collapses onto the surface of the 
secondary star during the close approach of the stars around 
periastron. This collapse would be initiated by the winds 
failing to attain a stable momentum balance at some point. 
Hydrodynamic simul ations of a co llapse of the WCR in 
the L Orionis system (|Pittardlll99i ) lead to a minimum in 
the fl ux and temperatur e of the X-ray emission at perias- 
tron (jPittard et al.ll2000l ). A similar event may occur in the 
rj Car system. In both the "collapse" and "discombobula- 
tion" scenarios mentioned above, it may take some consid- 
erable time to re-establish a bona-fide WCR. 



Table 3. The cfTect of increasing the primary star mass-loss rate 
on the prcshock wind speeds {vi and D2), the effective momentum 
ratio along the line-of-centres irn^^), the distance of the stag- 
nation point from the centre of the companion star (r'2), and 
the cooling parameter of the secondary wind (x2)- v\ and V2 
are calculated using /3-velocity laws (i)ooi = 500 kms"""^, v^2 = 
3000 kms-\ /3i = 4, and /32 = 1). M2 = 1.4 X lO-^Moyr"! 
in all models (as determined in § 14.2. Il l and the assumed stellar 
separation is dsep = 359 R© (corresponding to periastron in our 
standard model with e = 0.9). The radius of the gravitationally 
bound core of the primary star is taken to be 100 Rq. ?7ioc and 
X2 are evaluated using the preshock, rather than the terminal, 
wind speeds along the line-of-centres. The orbital aberration of 
the WCR is neglected for simplicity. 



Ml 


Vl 


"2 


Vloc 


r2 


X2 


(Msyr-i) 


(kms-i) 


(kms-i) 
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5.0 X 10"* 


127 
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0.58 
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2.0 X IQ-^ 


164 


2295 


0.10 


0.24 


1.2 


3.0 X 10-3 


173 


2130 


0.063 


0.20 


0.71 



5.1 Effects of wind acceleration and variable mass 
loss 

The models in §|4]have been made with the assumption that 
the winds collide at terminal velocity. In reality it is likely 
that one, or both, of the winds may still be accelerating when 
they collide. To gauge how this affects the previous results, 
the position of the stagnation point has been computed us- 
ing /3-velocity laws for both winds, with P — 4 and 1 for the 
primary and secondary winds respectively. Compared to the 
terminal winds case, ram pressure balance occurs closer to 
the primary star at all orbital phases (since its wind accel- 
erates more slowly). For instance, at = 1.0 the preshock 
primary wind only reaches a velocity of vi — 125 kms~^ and 
ram pressure balance occurs at a distance of 0.44 dsep from 
the companion star, compared to 0.30 dsep when terminal ve- 
locity winds are assumed. Cooling parameters for both winds 
were then calculated. As expected, the primary's wind was 
found to be highly radiative throughout the orbit (x ^ !)• 
Interestingly, the cooling parameter of the shocked compan- 
ion wind at periastron is at roughly the point that cooling 
starts to become non-negligible (x 3 — 5 see Fig. I16|l . 

Cooling of the shocked companion's wind becomes more 
important if there is a period of enhanced mass-loss from the 
primary wind (see Table [3| since the WCR is pushed deeper 
into the secondary wind resulting in lower preshock veloc- 
ities and higher densities (see Fig. I17|l . Note the dramatic 
change in the width of the postshock secondary wind (i.e. the 
distance between the shock and the contact discontinuity) 
on the line-of-centres as X2 becomes smaller. The intrinsic 
X-ray emission from these models is shown in Fig. 1181 where 
two trends are noticeable. First, the total flux increases with 
the primary mass-loss rate since a greater portion of the 
secondary wind is shocked, and second, the emission also 
becomes softer, since the preshock secondary wind speed 
declines. It is clear that the net effect of increasing Mi is 
to increase the intrinsic X-ray emission, though this would 
only be a temporary effect if any increeise in Mi was short- 
lived. Any increase in Mi cannot last too long as otherwise 
it would lead to an observed increase in the X-ray flux as the 
lines-of-sight to the emitting plasma exit the primary wind. 
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Figure 18. 1-10 keV spectra showing the variation in the intrin- 
sic emission calculated from a 2D hydrodynamic simulation as 
the mass-loss rate of the primary star increases. The spectra are 
labelled by their respective values of X2, as noted in Table [3] 



5.2 Radiative braking/inhibition 

The radiation flelds of both stars may affect the position and 
nature of the WCR around periastron passage. In this sec- 
tion we consider two possible effects: radiative braking and 
radiative inhibition. The former refers to the deceleration of 
a stellar wind by the radiative flux from the opposing star 
prior to reaching the WCR, whereas the latter refers to the 
reduction of the net rate of acceleration of a stellar wind 
due to the opposing radiation fleld. With this in mind it is 
useful to note that inhibition occurs close to the star driv- 
ing the wind, i.e. in the acceleration region of the wind, and 
braking occurs close to the WCR. To determine whether ra- 
diative braking is e ffective for either star, the equations of 
iGavlev et al.l (Il997l 'l have been evaluat ed. To perform these 
calculations the line driving theory of ICastor et al.l (119751 ) 
(hereafter CAK) was used to find the values of the param- 
eters k and a required to produce the desired terminal ve- 
locities and mass-loss rates for each star (Table [Jjl . While 
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Figure 17. 2D hydrodynamic simulations of the WCR at a separation corresponding to periastron for e = 0.9. The plots show the effect 
of increasing the primary mass-loss rate (Mi = 5.0 X 10~* (left panel) and 3 X 10~'^ (right panel) M0yr~^ ). In each simulation the 
wind speeds are set according to the on-axis preshock speed determined from a simple consideration of momentum balance (see Table |3] 
where the effective wind momentum ratio along the line-of-centres is also given). As the WCR is pushed into the acceleration zone of 
the secondary wind the shocked secondary wind becomes increasingly radiative (x2 = 3.6 in the left panel and 0.71 in the right panel) 
and the WCR becomes increasingly unstable. No gravitational or radiative driving forces are included in these calculations. 



Table 4. The parameters used to calculate the line driving of 
the stellar winds, and the effect of radiative braking/inhibition. 
Tcs is the temperature at the surface of the gravitationally bound 
core of the star, which for the pri mary star is t aken to be at a 
radius of 100 Rq. k and a are the lCastor et al" line driv- 

ing parameters, where subscripts 1 and 2 are used to define the 
coupling (C) between the winds and the radiation fields of the 
primary and secondary stars respectively. 





Primary 


Secondary 


TcB (K) 


25,800 


30,000 


L. (106Lq) 


4 


0.3 


k 


0.351 


0.443 


a 


0.525 


0.712 


M(M0yr-l) 


5.0 X 10"* 


1.4 X 10~5 


Voa (kms~^ ) 


500 


3000 



there is some uncertainty to the appropriate values of k 
and a these values provide a reasonable starting point for 
examining the dynamical interaction between the radiation 
fields and the winds in this system. We define coupling Ci 
as the the coupling between the primary star radiation and 
stellar wind, which is described by k\ and ai, and respec- 
tively coupling C2 constitutes the equivalent for the com- 
panion star. An un certainty in the use of the equations in 
iGavlev et all (|l997l ) is that it is not clear which values of k 
and a should be used when considering the braking of one 
star's wind by the radiation field of the other (i.e. should 
the primary star's k and a values be used when con sidering 
braking by the secondary star's radiation ? see also iPittardI 
1 19981 : [st -Louis et al.ll2005l . for further discussion). Therefore 



we have investigated a number of different possibilities. We 
find that the primary star's radiation does not brake the 
companion's wind, irrespective of which CAK parameters 
from Table |4] are used. When the close proximity of the 
WCR to the companion star is considered this is not such a 
surprising result as the ability of sudden radiative braking 
to halt a wind is expected to be effective when the incom- 
ing wind is close to impacting the surface of a star, which is 
not the case for the companion star's wind travelling towards 
the primary star. However, the opposite scenario of the com- 
panion's radiation field braking the incoming primary star's 
wind, does occur. Fig. [19] shows that in this respect alone 
at an orbital eccentricity of e = 0.90, braking is significant 
when the coupling is with either star (i.e. coupling Ci or C2). 
However, unless the coupling Ci is used, the braking does 
not occur before the primary wind has reached the WCR. 
The location of sharp braking moves closer towards the sec- 
ondary star as e increases, meaning that a normal ram-ram 
balance is increasingly likely. It is unclear whether the high 
density of the primary's win d will in reality preven t effective 
braking, as the argument of iKashi fc Sokeil (|2008h does not 
consider the ram pressure of the companion wind acting on 
the dense blobs of primary wind and the subsequent rate of 
ablation which may destroy the blobs before they reach the 
companion star's surface. A more in-depth examination of 
this mechanism is beyond the scope of this work, but will 
be considered in future models. Ou r current calculations im- 
prove upon the approach taken bv lKashi fc Sokeil (|2008l ) by 
determining the full range of coupling parameters. 

As mentioned above, the radiation fields of the stars 
can also play a prominent role in altering the mass-loss rates 
and terminal velocities of the stellar winds by inhibiting the 
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Figure 19. Radiativcly braked velocity profile of the primary 
star's wind as it approaches the companion star for various orbital 
eccentricities and coupling parameters. Radiative braking causes 
the pronounced sudden reduction in the velocity of the incoming 
wind (from the right in this plot). Orbital motion is not included. 
The parameters used in the line driving calculations are noted in 
Tabled 



acceleration of each other's wind (jStevens fc Pollo"c^ll994 ). 
When applied to t) Car we find that the radiation field of 
the primary star significantly reduces the terminal velocity 
of the companion wind (Fig. [20} to ~ 1520 kms^^ (assum- 
ing coupling C2). A further reduction occurs if the gravita- 
tional influence of the primary star on the companion star 
is considered. In this case tidal deformation results in the 
radius of the secondary star towards the primary increas- 
ing by ~ IR0. This in turn reduces the region available for 
wind acceleration, and the terminal velocity and mass-loss 
rate of the companion's wind drop down to ~ 1420 kms^^ 
and 1.25 x 10~^ Mq yr~^ respectively. When calculations are 
performed with coupling Ci, the companion's wind speed is 
significantly reduced at all orbital phases and the preshock 
velocity decreases such that a stable momentum balance 
cannot be attained at periastron (Fig. I21|) . 

These alterations to the wind parameters will have dras- 
tic implications for the observed e mission and may exp lain 
the wind disturbance discussed bv lMartin et al.l (|2006l ). To 
examine this scenario we have performed further hydrody- 
namic simulations where the wind speeds of the stars are 
vi_ — 360 kms^^ and V2 = 1420 kms^'^ (i.e. accounting 
for the slow acceleration of the primary wind, and the ra- 
diative inhibition of the secondary wind with coupling C2). 
The shocked companion wind reaches peak temperatures of 
< 10^ K in this case, compared to ~ 10* K if colliding at 
3000 kms~^. The intrinsic spectra calculated from these hy- 
drodynamic simulations are shown in Fig. [22l There is a 
dramatic reduction in the emission above 2 keV when the 
preshock velocity of the companion is reduced to the value 
implied by the radiative inhibition calculations. This pro- 
vides some explanation for the over-estimation of the hard 
band flux in the dynamic model in § [4] where terminal ve- 
locity winds are assumed (see Figs. ll4l and ll5|l . Models with 
reduced preshock companion wind velocities yield much im- 
proved fits to the XMM-Newton spectra in Fig. 1151 For in- 
stance, ai (j) = 0.924 a better fit can be attained when a 
preshock velocity of 2200 kms~^ is used, which indicates 



that radiative inhibition is important at this phase. Also, 
the excess in high energy flux seen in the 'standard' model 
spectra at (p = 1.009, 1.015, 1.018, and 1.023 can be rectified 
when reduced preshock companion wind velocities are used 
(a noteable failing, however, is the discrepency which still ex- 
ists between the model and the data at = 0.988 and 0.990). 
Interestingly, the required velocities lie within the range of 
predicted preshock velocities caused by radiative inhibition 
(Fig. I21|l . Considering the recovery of the companion wind 
as the stars recede, we would expect a preshock velocity of 
V2 — 2500 kms~^ at ~ 1.13, which is significan t ly hig her 
than the ~ 1100 kms"^ detected bv llping et"aLl (|2005l ) at 
this phase. However, this comparison does not account for 
projected velocity vectors and/or the preshock velocity of 
the downstream gas. We note that although there is some un- 
certainty in the coupling parameters used, the calculations 
performed provide a useful guide to the range of inhibited 
velocities. 

It is also useful to consider the degree of radiative in- 
hibition near apastron. PC02 found a best-fit velocity of 
3000 kms^^, and this value has been adopted in our 'stan- 
dard' model. Reasonable agreement with this value is at- 
tained when using coupling C2, where we find the termi- 
nal wind speed is reduced to Vcx,2 = 2916 kms~^. In con- 
trast, with coupling Ci, the companion's wind speed is more 
than halved giving Voc2 = 1321 kms~^. Hence, our calcu- 
lations support a coupling which is inclined towards the 
wind rather than the radiation field in qu estion. At first 
sight this is in contrast to the findings of ISt-Louis et al.l 
(2005|), where the wind- wind collision in Sanduleak 1 was 
investigated, although their conclusion may depend on the 
assumed mass-loss rates. It should be noted that in reality 
the values of k and a are spatially and phase dependent 
(e.g. IPuIs et al.|[2006l : iFuUerton et aLlliooil ). yet to make the 
problem tractable we needed to make some simplifying ap- 
proximations. 

Clearly, radiative inhibition is an important mechanism 
in the rj Car system, and its effect must be properly consid- 
ered in future models. 



5.3 Collapse of the wind collision region onto the 
secondary star 

Examining the ram pressure balance between the reduced 
velocity winds discussed in the previous section reveals that, 
when e = 0.9, the equilibrium is stable when the compan- 
ion's k and a (i.e. coupling C2) are used in the radiative 
inhibition calculations (Fig. I23|l . An increased orbital eccen- 
tricity reduces the ram pressure of the companion's wind 
at periastron, due to the increasing influence of the pri- 
mary's radiation fleld as the stellar separation is reduced. 
At e = 0.95 there is no stable balance between the winds 
at periastron, which is indicated by the lack of an intersec- 
tion between the ram pressure curves for the winds with 
opposing gradients. If coupling Ci is used to perform these 
calculations there is no stable balance point at periastron 
even when e = 0.90. Clumping within the winds, a tempo- 
rary increase in the primary mass-loss rate, and the gravi- 
tational influence of the companion star may further reduce 
the likelihood of a stable balance between the winds. Refer- 
ring back to Fig. [19] it is clear that in the scenarios where 
a collapse is predicted, subsequent radiative braking of the 
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Figure 20. The speed of the companion wind towards the pri- 
mary star (located at a distance of 359 which corresponds 
to the separation of the stars at periastron for e = 0.9) assum- 
ing that the coupUng of the primary star's radiation field to the 
companion star's wind is described by coupling C2 (i-e. the same 
values of the CAK parameters as required for the companion star 
to drive its own wind). Clearly, the acceleration of the wind is 
strongly inhibited by the radiation field of the primary star. 
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Figure 21. The preshock velocity of the companion's wind as a 
function of orbital phase for e = 0.9. Three scenarios are plotted: 
no radiative inhibition (no Rl), Rl with coupling C2 (see ii 15. 2t , 
and Rl with coupling Ci. The preshock velocity was calculated 
after determining the location of the stagnation point between 
the winds. The lower curve does not have a stable balance point 
at periastron, and the companion cannot launch a wind towards 
the primary star. Orbital motion is not included in this plot. 
The parameters used in the line driving calculations are noted 
in Table |4] The squares represent the preshock velocities used to 
attain improved fits to the XMM-Newton spectra in Fig. 1151 



primary star's wind before it reaches the companion star's 
surface should also occur. This is of course a very compli- 
cated picture, and in reality it is difficult to predict the sum 
of all of the different mechanisms described above without 
a single model which self-consistently includes all of the ap- 
propriate physics which we have shown is important. This is 
beyond the scope of the current work, but will be addressed 
in future models. 

Is there any observational evidence for a collapse of 
the WCR onto the companion star? Let us first consider 
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Figure 22. Intrinsic spectra calculated from hydrodynamic simu- 
lations discussed in ii l5.2l f corresponding to periastron for e = 0.9), 
where the winds are assumed to collide at terminal velocities 
(solid) or at much reduced velocities (dashed). 



the X-ray emission. While radiative inhibition can dramati- 
cally reduce the hard X-ray emission around periastron (see 
Fig. I22|l , it cannot by itself explain the long duration of the 
X-ray minimum and its slow recovery, since the stars sepa- 
rate very quickly after periastron. Furthermore, the preshock 
wind speed of the companion inferred from the fits to the 
XMM-Newton data remains roughly constant during the X- 
ray minimum {(f) ~ 1.0 — 1.03), whereas the inhibition cal- 
culations indicate that it should rapidly rise (Fig. 1211) . 

However, if the WCR collapses onto the surface of the 
companion star it is unclear how quickly the companion 
wind will be able to re-establish itself. It is entirely possible 
that once the primary wind collides directly with the surface 
of the companion star its continued ram pressure will pre- 
vent the companion star from developing a wind towards the 
primary star until long after periastron has passed (see the 
hydr odynamic simulations of such an eventuality in IPittardI 
Il998l ). Indeed, the asymmetry of the X-ray minimum about 
periastron passage indicates that the WCR in 77 Car exhibits 
hysterisis in this sense. 

To further examine this possibility a mock collapse of 
the WCR has been incorporated into the dynamic model in 
§ m The collapse is assumed to begin at = 1.0 and to 
end at = 1.03. The emission from the region of the WCR 
which is affected by the collapse is set to zero, while shocked 
gas further downstream (and upstream after the recovery) 
emits as normal. The simulated observed emission following 
this addition to the model is displayed in Figs. [T4l 1151 and 
1241 The reduction in the hard band emission relative to the 
previous models is clearly evident in the spectra shown in 
Fig. 1151 and a much better match to the observations during 
the minimum compared to the 'standard' model is obtained. 
The fit to the lightcurve minimum is also considerably im- 
proved (Fig. I24|) . While there are remaining differences in 
the hardness ratio (Fig. I14p . the similarity of its shape dur- 
ing the X-ray minimum is intriguing (specifically the slower 
change from softer to harder emission in the second half of 
the minimum). Referring back to our conclusion in 14.1.21 
that lines-of-sight which favour the companion star in front 
at periastron are excluded still remains true in the collapse 
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Figure 23. The ram pressure (top), velocity (middle), and den- 
sity (bottom) of the radiatively driven winds as a function of 
distance from the companion star and orbital eccentricity. The 
separation of the stars corresponds to periastron in each instance. 
In each plot the companion's wind (thin lines) starts at 20 Rq. 
Radiative inhibition of the companion's wind by the primary's ra- 
diation field is included, where coupling C2 is used. The top panel 
shows whether a stable momentum balance can be achieved (in- 
dicated when an intersection occurs with opposing gradients, sec 
Fig. 2 of IStevens et al.lll992l) . The braking radius (Fig. [I9j for 
each scenario is indicated by the marks at the top of the upper- 
most plot. No stable balance exists for e > 0.95 unless braking is 
important. The parameters used in the line driving calculations 
are noted in Table |4] 
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Figure 24. The 2-10 keV X-ray lightcurvc showing the effect of a 
collapse of the WCR onto the companion star between phases 1.0 
and 1.03. During this time the only emission arises from further 
downstream. The other parameters of the model are the same as 
those for the r] = 0.18 lightcurve shown in Fig. 1131 The crosses 
represent the integrated fluxes from the individual models with 
lower companion wind velocities discussed in ij 15.21 



scenario as the pre/post-minimum fiux level provides a poor 
match to the observed lightcurve. 

Therefore, these results provide significant support for 
a collapse of the WCR onto the companion star around pe- 
riastron. 

During the collapse of the WCR the com panion 
star may accrete some of the primary star wind (jSokeil 
Taking the mass accre tion rate to be M ace — 
■n:pi{r)Rl^^Vi{r) (e.g. Eq.(15) in lAkashi et al.ll2006l ), where 
pi(r) and ui(r) are the density and velocity profiles calcu- 
lated for the primary star wind and the accretion radiufl 
Rb.cc = R2 (Table [TJ. We find the total accreted mass dur- 
ing the collapse to be ^ 7 x 10~* Mq , which is approximately 
one to two ord e rs of magnitude lower than values noted by 
iKashi fc Sokeil (|2008l ). 



6 CONCLUSIONS 

We have constructed the first 3D model of the colliding 
winds in rj Car with spatially extended, energy dependent, 
X-ray emission. Coriolis forces skew the WCR and wrap 
the far downstream regions around the stars. Attenuation of 
the X-rays through the unshocked winds of both stars, and 
through the cold, dense layer of postshock primary wind is 
considered. The simulations place constraints on the orbital 
orientation and nature of -q Car through comparison with 
observations from the RXTE and XMM-Newton satellites. 
The main findings from this work can be summarized thus; 

• 3D effects, specifically the skew to the WCR due to or- 
bital motion, lead to an asymmetry of the X-ray luminosity 



^ In lAkashi et al.l 1I2OO6I ) the accretion radius is determined under 
the assumption that the companion star is not driving a wind in 
any direction. In contrast, we assume that the companion star 
drives a wind on the side facing away from the primary star, 
hence accretion only occurs for gas striking the near-side of the 
companion star and the accretion area becomes ~ 7ri?| 
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either side of the minimum. Circumstellar absorption is usu- 
ally dominated by the un-shocked stellar winds. For a very 
brief period [Scf) ^ 0.01) the absorption may be dominated 
by material in the cold postshock primary wind. 

• Inclination angl es similar to the polar axis angle of the 
Homunculus nebula (|Smithll2006l ) are favoured. 

• Orientations with the secondary star in front of the pri- 
mary around periastron passage are supported. The model 
results favour a line-of-sight angled at — 30° to the semi- 
major axis in the prog rade direction, in ag reement with 
iNielsen et"aLl (|2007h and lOkazaki et al.1 (|2008i l. Orientations 
with the secondary star in front of the primary or at quadra- 
ture during periastron passage are excluded. 

• With terminal wind velocities fixed at Vooi = 
500kms~^ and Vaa2 = 3000 km s~^ for rjCax and the com- 
panion star respectively (our 'standard' model), the best 
fit to the 2-10 keV RXTE lightcurve was obtained with 
mass-loss rates of Mi = 4.7 x 10"''MQyr~^ and Ah = 
1.4x 10~^ Mq yr~^ . These values are similar to an earlier de- 
termination by PC02 though the mass-loss rates are slightly 
higher and the wind momentum ratio, r] is slightly lower. 
This difference reflects changes in the adopted distance and 
interstellar absorption to the system, and to difficulties in 
separating other, more spatially extended, components from 
the large RXTE beam. It may also reflect the lack of mixing 
at the CD in the dynamic model used in this work, compared 
to the hydrodynamical simulations of PC02. 

• Significant discrepancies between the observed and 
'standard' model lightcurves and spectra exist through the 
X-ray minimum in. We conclude that a wind eclipse is not 
the sole cau se of the minimum, in contrast to the recent con- 
clusions of lOkazaki et al.l (|2008l 'l who modelled the X-ray 
emission as a mono-energetic point source and considered 
only the lightcurve. 

• An examination of the importance of radiative inhibi- 
tion has revealed that around periastron the preshock ve- 
locity of the companion wind is likely to be significantly 
reduced. The secondary wind is then shocked to much lower 
temperatures and cools rapidly. As a result the hard X-ray 
flux is efficiently quenched. All future models must properly 
account for the effect of radiative inhibition. 

• While the hard X-ray flux can be significantly reduced 
during periastron passage, the long duration of the mini- 
mum after periastron when the stars are rapidly separating 
hints that something very extraordinary happens. We be- 
lieve that the most likely explanation is the collapse of the 
WCR onto the companion star, perhaps helped by a short 
increase in the mass-loss rate of the primary wind. Incorpo- 
rating a mock collapse into the model between phases 1.0 
and 1.03 yields significantly better agreement between the 
model and data. 



The dynamical model described in this paper provides 
insight into the role of 3D effects in highly eccentric bina- 
ries. Its failings to produce synthetic X-ray emission which 
is a good match to the observed data have highlighted the 
importance of the stellar radiation fields on the dynamics 
of the winds and the structure and cooling of their collision 
region. To advance this work further requires 3D hydrody- 
namical models which incorporate cooling, gravity and the 
driving of the stellar winds. 
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